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l have worked on innumerable problems that you would call humble, but which I enjoyed and felt very good 
about because I sometimes could partially succeed . . , Richard P. Feynman, letter to Koichi Mano 

INTRODUCTION 

Harold Morton introduced a talk by saying that when you wind up an old professor, he tends to talk for a 
microcentury. I will attempt to keep my comments to that canonical time span. Having failed to find some 
unifying theme for this talk, I decided to just ramble through my career with a focus on the algorithms, 
spacecraft, and people I’ve had the privilege and pleasure to work with. The algorithms, and certainly the 
spacecraft, are not all mine. The people are some of those whose ideas that have most influenced and 
inspired my career. The organization of the paper is largely chronological, but I do not hesitate to jump 
forward or backward in time when the material demands it. The coverage is broad but necessarily shallow; 
the interested reader can find more detail in the references. 

PREHISTORY 

My father was an engineer, and engineering was a natural career path for my older brother and me. I much 
preferred arithmetic to other courses in school, because arithmetic problems have right and wrong answers. 
It was not until late in my undergraduate education that the realization dawned on me that engineering, and 
research in general, required dealing with uncertainty and ambiguity. However, the instinct to search out 
exactly solvable problems has never left me. 

When I was an adolescent, I found an article on Calculus in an 
encyclopedia. It contained the well-known problem illustrated 
at the right. A tinsmith is to cut corners out of a square of 
metal and fold up the edges to make an open box of maximum 
volume. The optimum cut is given by the solution of 

0 = d\x(a- 2x) I 2 \/dx =(a-6x)(a-2x) . (1) 

This is a humble problem indeed, but I found it incredibly 
exciting that a “practical” problem could be solved by simple 
algebra. One could say that this marks the beginning of my 
interest in optimization problems, and I’ve never really lost 
my sense of wonder at the power of mathematics. Figure 1 . First Humble Problem 

I followed my brother to Cornell University, but chose to study Engineering Physics rather than Chemical 
Engineering as he did, mostly because it was reputed to be the most challenging major. The Cornell 

Mathematics Department offered a five-semester course in applied mathematics for Engineering Physics 
majors, taken after the usual three-semester calculus sequence. This was an outstanding sequence taught in 
small sections not by teaching assistants but by faculty members, including the brilliant applied 
mathematician Mark Kac, who taught differentia] equations to sophomores. In my ninth semester of the 
five-year Engineering Physics program, I elected to take a course in probability theory taught by Jacob 
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Wolfowitz, whose son currently heads the World Bank. An undergraduate course in quantum mechanics 
seduced me away from Engineering Physics to pure physics, and I pursued graduate studies in theoretical 
physics at the University of California at Berkeley. I spent a wonderful five years there, with future Nobel 
laureates as teachers and one colleague. In retrospect, though, the only practical things I took away from 
Berkeley were a familiarity with group theory, the only useful mathematics I hadn’t learned at Cornell, and 
the experience of writing a short FORTRAN program for my dissertation. 

After Berkeley, I was' awarded a postdoctoral fellowship at the University of Maryland. My principal 
achievements there were to meet Malcolm Shuster, who was a graduate student in physics at the time, and 
to reinforce my conviction that theoretical physics was not the field in which I would make my mark. 
Believing that undergraduate teaching would be more to my taste, I declined a second year of my 
fellowship to accept a position as Assistant Professor of Physics at Williams College. 

It’s been said many times that you never really learn a subject until you teach it, and my years at Williams 
certainly bore this out. Aside from deepening my knowledge of physics, teaching at an excellent small 
liberal arts college also fostered interactions with mathematicians, other scientists, and humanities faculty. I 
greatly enjoyed this life, perhaps too much, for I did not earn tenure. The spring of 1974 found me at the 
Washington meeting of the American Physical Society, searching with many others for a congenial 
teaching position. As luck would have it, the Computer Sciences Corporation (CSC) was at the meeting to 
recruit scientists with experience in computer programming to write software for NASA. My resume 
included a course on Use of the Computer in the Physical Sciences that I had taught on an IBM 1 130 at 
Williams. This helped earn me an interview with CSC, and I was offered a position one month before I 
would have become unemployed. 

COMPUTER SCIENCES CORPORATION (CSC) 

Quaternions and Kalman Filters 

CSC offered me a position in either the Attitude Systems Operation or the Orbit Systems Operation. I chose 
the former, mostly because I felt a rapport with the two physicists in that group who had interviewed me: 
Jerry Lerner, who had been Malcolm Shuster’s roommate in graduate school, and Dale Headrick. Dale was 
the first to arouse my interest in quaternions and Kalman filters. I had never heard of quaternions, since 
they did not appear in any of the five classical mechanics texts I had learned or taught from. Goldstein did 
not introduce the quaternion components as the real and imaginary parts of the Cayley-Klein parameters 
until the second edition of his classic text [1], I had heard of Kalman filters, though. The math department 
at Williams offered annual informal discussion groups of selected topics of interdisciplinary interest. 
Ironically, the topic chosen for my last year at Williams was Kalman Filtering; but I was too busy looking 
for a job to take advantage of that opportunity. 

My first supervisor at CSC was Mike Plett, who said that he didn’t want to see any clever programming 
from me. I thought that was very peculiar advice until I wrote a clever program that inadvertently assigned 
a new value to the constant 1. After that, I concentrated on good, solid programming, following the 
example of Myron Shear, who was described by Dale Headrick as the best FORTRAN programmer he’d 
ever known. 

Attitude/Orbit Analysis 

My first task at CSC was under the direction of Mel (Tom) Velez at Goddard Space Flight Center (GSFC). 
He was pursuing simultaneous orbit/attitude estimation of the spin-stabilized Synchronous Meteorological 
Satellite (SMS), since he realized that Earth horizon sensor data contains orbit information. This led to 
some interesting work, but the vastly different time scales of the orbit and attitude dynamics impeded 
progress. The orbital equations of motion could be integrated with time steps of tens of minutes, while time 
steps of fractions of a spin period were required by the rigid-body attitude dynamics. This led Velez to 
pursue variation-of-parameter (VOP) approaches to rigid-body dynamics. Harold Morton, John Junkins, 
and their students had done some interesting work in this area, but I tried a different approach [2, 3], 
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The usual approach was to integrate the equation obeyed by the inertial-to-body attitude quaternion 


dq 

dt 


®q 


( 2 ) 


and the equation of motion for either the inertial components of the angular momentum 

dL I /dt = N I 


or for the body-frame components 


dL B /dt = N 3 — to x L 5 , 


(3a) 

(3b) 


where to denotes the angular velocity and N the external torque. This approach has one redundant equation, 
which is constrained by the unity norm of the attitude quaternion. The problem with this formulation is that 
all the quaternion components are “fast” variables unsuitable for a VOP approach. My idea was to use both 
of Eqs. (3a) and (3b) and not Eq. (2). There again is one constraint, since the norm L of the angular 
momentum is the same in both frames. One more variable is needed to complete the parameterization, 
which I chose to be an angle variable %, obeying the “fast” equation 


d X /dt = L(L 2 + L^ •L / )~ 1 [(L 8 + L 7 ) ■©+ r 2 (L B xL / )-(N b + N / )] . (4) 


Unfortunately, the VOP approaches to attitude dynamics never bore fruit, and these equations (in a much 
cruder form) were put away in a drawer. Many years later, they found a home in the special issue on 
Attitude Representations of The Journal of the Astronautical Sciences edited by John Junkins and Malcolm 
Shuster [4], These variables have never found widespread use, but they have recently been employed for 
attitude estimation [5], 

Gene Lefferts 


After about a year of this attitude/orbit analysis, Velez moved on to other matters, and my task fell under 
the supervision of Eugene J. Lefferts. This was the most felicitous transition of my career. Gene Lefferts 
became my most important mentor, and I owe my success from that point onward to his influence. He was 
also directly responsible for two of my future career moves. When I first began working with Gene, one of 
my CSC associates advised me to talk to him about fishing, because that was his interest outside of work. I 
quickly found that my complete lack of interest in fishing was no drawback, because Gene was interested 
in everything: politics, music, art, and culture in general, as well as fishing. He advised an Explorer post 
that explored fishing, caving, avant garde theater, and gamelan music, among many other things. My 
daughter joined Gene’s post as an adolescent, and he was an important influence on her as well. 

Gene actively promoted the use of quaternions, Kalman filters, and dynamic time-domain simulation; and it 
is owing to his influence that these have been major themes in my career. One of Gene’s insights was that 
optimal estimation, in the form of the Kalman filter, had been very important, but that optimal control had 
never demonstrated any superiority to classical control techniques. This insight will sound very strange to 
the astrodynamics community, where optimal maneuver planning has been invaluable in space mission 
design. I believe that Gene was correct in the context of attitude control, where the control effort can be 
generated with little or no expenditure of fuel. 

Gene also helped me develop my intuitive feel by means of humble toy problems. One example is the 
simple problem of finding the steady-state covariance of a one-dimensional state with constant dynamics 
and measurements processed at a time interval A t . The pre- and post-update covariances are given by 

P(-) = P(+) + QAt (5a) 

P(+) = [I-KH]P{~) with K = P(r)H T [HP(-)H T + RY 1 (5b) 
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We can assume H = 1 and R = a 2 , and the multiplication of lxl matrices is commutative, so these 
equations are easily solved to give 


^(QAtf + ^QAt + QAt^ and P{+) = \ yjiQAt) 2 + 4a 2 n QAt - QAt 
Filtering is only really useful when QAt <s. R , in which case 


P( -)=2 


P(-)~P(+)~(J„jQAt. 


( 6 ) 

(7) 


These equations are obviously useful in giving the expected estimation accuracies given specified levels of 
measurement noise R and process noise Q. Gene pointed out to me the less obvious use of these equations 
to determine the level of dynamic measurement fidelity, as represented by Q, needed to reduce estimation 
errors to a required level given specified measurements. It’s interesting to note that Eq. (7) is more easily 
derived as the solution of the differential equation arising from the continuous-measurement limit of the 
Kalman filter [6], which gives the steady-state covariance as the solution of 

dP/dt = Q-PH T (RAt)~ ] HP = 0. (8) 

I have gained many very useful insights by looking at the steady-state behavior of simple Kalman filters. 

Working with Gene, I developed a general quaternion-based simulation that accommodated a variety of 
sensor models and control strategies. We used this in an attempt to use dynamics modeling to improve the 
attitude estimation of the Nimbus spacecraft [7], This was of interest to Harry Stallings, an engineer in 
Henry Hoffman’s Guidance and Control Branch at Goddard who had been a colleague of Gene’s atthe 
Martin Aircraft Corporation in Baltimore. The Nimbus analysis was inconclusive, but Gene and I had 
established a connection that was to prove important when Goddard began development of the Solar 
Maximum Mission (SMM) spacecraft. 


Solar Maximum Mission (SMM) 


SMM was the first of the Multimission Modular Spacecraft (MMS) series, and was later to be the first 
spacecraft to be serviced in-orbit by astronauts. The MMS had three modules: a power module, a 
communication and data handling module, and the Modular Attitude Control System (MACS). The first 
analysis we performed for SMM was a dynamics simulation of the despin of the spacecraft after release 
from its launch vehicle, using magnetic torquers. We showed that despin from worst-case conditions could 
require ten hours, although the initial analysis had given a time of one hour. It turned out that the earlier 
analysis had misplaced a decimal in one of 
the conversions of magnetic moment from 
Ampere-ft 2 to Ampere-m 2 , magnetic field 
strength from gammas to nT, and the 
resulting torque from oz-in to Nm. This 
was an early lesson in the advantage of 
using SI units consistently, as we did in our 
dynamic simulator, with a penalty much 
smaller than that paid in 1999 by the Mars 
Climate Orbiter. Our analysis also led to 
providing SMM with an additional 
magnetic torquer bar external to the 
MACS. I was delighted that my 
contribution had actually led to a change in 
hardware launched into space. 



Figure 2. SMM with Attached Astronaut 
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Jim Murrell, of General Electric Space had developed a multimission quaternion Kalman filter for the 
MACS being built by GE [8], This Kalman filter was later used on Landsats D and E and on the Upper 
Atmosphere Research Satellite (UARS); but Arun Guha, Henry Hoffman’s Lead ACS Engineer for SMM, 
chose not to use it for SMM. Instead, he developed a decoupled-axis Kalman filter based on the filter 
employed on GSFC’s International Ultraviolet Explorer (IUE) spacecraft. With our credibility established, 
Gene and I were given the task of developing a FORTRAN prototype of the ACS software for SMM and 
validating it in our dynamic simulator. This was the most exciting project I had engaged in thus far, and I 
attacked it with great enthusiasm. After the prototype had been developed and testing was underway, Henry 
Hoffman decided that his branch should have a peer review of our work. The review panel included Jim 
Donohue, an outstanding classical controls engineer, also from the Martin Aircraft Corporation, who 
mentored an entire generation of engineers in Hoffman’s branch. He had also one time played on a semipro 
baseball team with the future Hall of Fame baseball player A1 Kaline, which impressed me very much. 
During my presentation, Donohue suggested that we apply a lead-lag filter to the digital sun sensor data. I 
had no idea what he was talking about but said that we’d look into it. When I later confessed my ignorance 
to Guha, he advised me to read the Classical Controls book in the Schaum’s Outline series so I could speak 
the language [9]. I took Arun’s advice, and later told him that I was taking a course in control theory. When 
he asked where, I said that it was the course that he was teaching me. I presented my SMM analysis at the 
same GN&C conference where Murrell presented his MACS Kalman filter [10], 

SMM was launched in 1980, after I had left CSC, and within a few months the fuses on the reaction wheels 
in the MACS began to blow, causing a loss of attitude control. Henry Hoffman and Tom Flatley, an 
extraordinary attitude dynamicist in Henry ’s branch, used SMM’s magnetic torquers to spin-stabilize the- 
spacecraft until astronauts could replace the entire MACS in 1984 (Figure 2). Henry liked nothing better 
than rescuing a satellite in distress, which led to his being called a satellite savior [11], 

Spacecraft Attitude Determination and Control 

During my stay at CSC, their contract for support of GSFC neared its expiration date. Roger Werking, who 
was responsible for the greater share of CSC’s support, tasked Jim Wertz to document the algorithms and 
software that CSC had developed under their contract. He intended this to be a Contractor Report that 
would be available for the use of another contractor that might outbid CSC for the continued support. Jim 
somehow talked Roger into producing a book rather than a Contractor Report, convinced CSC to contribute 
to its creation, and found a publisher. The resulting book has been an indispensable reference in spacecraft 
attitude determination and control and a monument to Jim’s initiative and perseverance [12]. 

Jim developed an overall outline for the book, and solicited the staff of the Attitude Systems Operation for 
proposals to write various sections. I offered to write the sections on attitude kinematics and dynamics, 
assuming that I only needed to paraphrase the corresponding sections of Goldstein [1]. You can imagine 
my chagrin when Mike Plett said, in reviewing my proposed outline of these sections, “I hope you can 
make this clearer than Goldstein.” I apparently succeeded, and found my small first taste of fame. Malcolm 
Shuster joined the CSC Attitude systems operation in 1977, too late to contribute significantly to 
Spacecraft Attitude Determination and Control. During his first year at CSC, he developed the QUEST 
algorithm, to which I had made some small contributions. This algorithm was announced to the world at the 
same GN&C Conference where Murrell’s and my papers appeared [13]. The history of the development of 
QUEST was the topic of Shuster’s Brouwer Lecture [14]. 

After the completion of Spacecraft Attitude Determination and Control, many of the authors began to leave 
CSC, owing partly to their employment opportunities increasing along with their visibility, and partly to 
their desire to avoid the uncertainties of employment by a government contractor. Ironically, CSC was to 
retain its supporting role at GSFC into the next millennium. I also became impatient to leave, preferably to 
enter the Civil Service. To my delight, an ad appeared in the April 23, 1978 Washington Post, seeking an 
Aerospace Engineer/Physicist/Mathematician to conduct “research in satellite orbit and attitude dynamics 



and the mathematical modeling of dynamical systems” at the Naval Research Laboratory. I applied with a 
cover letter stating, “I believe that my qualifications, as presented in the enclosed resume, are perfectly 
suited to this position.” My wife found this breathtakingly arrogant, but I thought it was a simple statement 
of the fact that the ad seemed to be written specifically for me. It turned out that Bernie Kaufman, who was 
seeking to fill this position, had worked with Gene Lefferts at Goddard and was checking out all the 
Goddard-related applicants with him. I later found that when Bernie told Gene that I had applied, Gene told 
him that he could stop looking further. This was the first of the two career moves I owe to Gene, 

NAVAL RESEARCH LABORATORY (NRL) 

Optimal Control 

Bernie Kaufman’s section comprised himself, Bob Dasenbrock, Bill Harr, and me. Bernie and Bob were 
the orbit dynamics experts, and Bill was an excellent programmer who showed me how to make my 
FORTRAN readable. Bernie is the best supervisor I ever had, supportive, encouraging, successful in 
bringing interesting work into our section, and insulating us from outside interference. 

Ironically, Terry Alfriend, who headed a different branch at NRL, posed the first problem that really 
engaged me. It was the rest-to-rest attitude maneuver about the symmetry axis of an axially symmetric 
spacecraft with one flexible mode. The equations of motion are 

I b 6=N + k\i/ and I a (8 + y/) = -ky / , (9a) 

with the boundary conditions 

8(± T/2) = ±9/2, 8(± T/2) = y/(± T/ 2) = y(± T/ 2) = 0 . (9b) 

where 0is the angular, position of the main body with moment-of-inertia I b , i//is the relative angular 
displacement of the appendage with moment-of-inertia I a , N is the torque applied to the main body, and k is 
the torsional coupling constant. I set up an optimal control problem with quadratic penalties on the control 
and the flexible mode. The problem has closed-form solutions with 6, t//and N all odd functions of time. 
When the penalty on the control was less than some critical value, the solutions involve hyperbolic trig 
functions. Replacing the quadratic penalty on the control with an absolute limit on its magnitude gives a 
more complicated solution with a singular control arc and only a slight reduction in maneuver time. 
Additionally removing the penalty on the amplitude rate of the flexible mode led to chattering arcs that 
defied analysis. At this point, I realized that, although the control effort became a more complicated 
function as I tried to improve the performance, Sand yr remained smooth functions of time. This led me to 
consider a simple non-optimal maneuver, choosing 8, \ff and N to be the lowest order polynomials that 
could satisfy the equations of motion and the boundary conditions. The maneuver time using this 
polynomial profile was only about 10% longer than the maneuver time of the optimal profile. This 
reinforced my agreement with Gene Lefferts’ belief that optimal control was less useful than optimal 
estimation. I presented this paper at the 1979 AAS Astrodynamics Conference, the first AAS conference I 
attended [15]. 

The Lefferts, Markley, Shuster (LMS) Paper 

Shortly after I moved from CSC to NRL, the analysis of steady-state Kalman filters regained my attention. 
The state equation for single-axis attitude estimation using a gyro in “model-replacement mode” is 



where 8 is the rotation angle, co s is the gyro output, b is the gyro bias, an n v and n„ are uncorrelated white 
noise processes. Assuming effectively continuous angle measurements, the steady-state covariance obeys 

dP/dt =FP+PF T + Q- PH T (RAt )~ 1 HP = 0 , (11) 
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with H = [1 0] , R = a 2 , and Q = diag([aj,a^]) . The solution, found without too much effort, is 

^=-<W,(A0 V2 , P ee = a n [{a 2 v -2P eb m 112 -o y a n {^ 2 , and P bb = o u {g 2 v -2pJ 2 ~o v o u . (12) 
The approximate form of P gg is basically the same as Eq. (7). 

The corresponding equations for discrete measurements lead to a quartic equation that is not so easy to 
analyze. Amazingly enough, Robert L. Farrenkopf found a factorization of this quartic into two quadratic 
equations, which led to an incredibly useful estimate of the accuracy of attitude estimation systems. His 
result was published in a note in the July-August 1978 issue of the Journal of Guidance and Control, right 
at the time I was moving from CSC to NRL [16]. 

I decided to try to generalize Farrenkopf’s result to the three-axis case, starting with the quaternion 
equivalent of Eq. (10). It soon became apparent that no closed-form steady-state solution would be 
forthcoming, but I persisted in developing the general equations. After much analysis, I arrived at a 
formulation that looked familiar; it was identical to the MACS filter developed by Murrell. Greatly 
disappointed, I sent my analysis off to Gene Lefferts at GSFC, asking if he thought it had any value at all. 
He said that Shuster, who had taken over my task after I left CSC, had developed a parallel analysis, and 
that he thought that our combined work was worth publishing. The resulting paper is the most cited paper 
that Malcolm Shuster or I ever wrote [17], The final version was written on weekends at Gene’s house. He 
provided the vector measurement model, the only real innovation in the paper, and also served as a buffer 
to prevent Malcolm and me from coming to blows over disagreements in presentation. 

Autonomous Navigation — • - - j 

My principal research activity at NRL was a covariance analysis of autonomous navigation using 
landmarks. This coupled attitude/orbit problem used some of the expertise I had developed with Velez, 
which had led to my being hired by NRL. The attitude and orbit errors turn out to be highly correlated, and 
it is not surprising that they tend to cancel when geolocating unknown objects on the Earth. My first paper 
on this subject was presented at the AAS Astrodynamics Conference at Lake Tahoe [18]. Connie 
Carrington, then a graduate student of John Junkins, presented a very nice paper on magnetic momentum 
unloading at the same conference. When I introduced myself and complimented her on her paper, she said 
that a paper Peter Camillo and I had written on our momentum unloading analysis for SMM had been very 
useful to her, which greatly pleased me [19], 

Further work on autonomous navigation augmented the landmark measurements with measurements of the 
line-of-sight (LOS) between two spacecraft in a formation. Remarkably, I found that the orbits of both 
spacecraft could be obtained using only measurements of the magnitude and inertial direction of the LOS 
between the spacecraft, except in a few highly symmetric situations [20]. This result was rediscovered 
independently by Mark Psiaki at Cornell [21]. 

As part of this covariance analysis, I developed a convenient form for the orbit state transition matrix 
expressed in Cartesian coordinates [22], Bob Dasenbrock later told me that this algorithm was used in the 
onboard software of at least one NRL spacecraft. 

Although the research atmosphere at NRL was very congenial, I missed the connection with real spacecraft 
that I had experienced in my work on SMM. In 1985, Gene Lefferts announced his retirement. Jerry Teles, 
his supervisor in the Flight Dynamics Analysis Branch, asked if there was a replacement who could supply 
the capabilities that he had provided to GSFC. He said that I could, and the position was offered to me. It 
was impossible for me to refuse, as hard as it was to leave Bernie Kaufman and NRL. Once again, I owed a 
career move to Gene Lefferts. 
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GODDARD SPACE FLIGHT CENTER 

Attitude Determination Error Analysis System (ADEAS) 

My first really interesting project after moving to GSFC was developing the Attitude Determination Error 
Analysis System (ADEAS). This was an attitude error covariance analysis program facilitating the analysis 
of the accuracies attainable with different sensor complements, different observation schedules, and 
different estimation methods, including mistuned estimators. It was developed for a mainframe computer, 
but has since been ported to the desktop computer environment and is still widely used at GSFC. One of 
my colleagues in this effort was Ed Seidewitz, a brilliant recent MIT graduate. Ed and I worked through the 
covariance analysis equations for both batch estimators and Kalman filters with solve-for and consider 
parameters. One of the innovations of ADEAS was its ability to assess the influence of process noise in a 
batch estimator. Mark Nicholson, the best FORTRAN programmer I have ever known, did most of the 
programming. He also uncovered errors in the equations Ed and I gave him to code. I believe that much has 
been lost in the increasing separation of the functions of computer scientist and control engineer. The 
results of Ed’s, Mark’s, and my seamless collaboration were documented in the best conference paper that I 
never published in an archival journal [23], One interesting discovery of our ADEAS work was that the 
attitude of Goddard’s Cosmic Background Explorer (COBE) spacecraft could be much more simply 
represented by 3-1-3 Euler angles than by any other representation of the attitude, weakening my aversion 
to this parameterization. 


Wahba’s Problem 

In 1965, Grace Wahba had posed the problem of finding the proper orthogonal matrix that minimized a 
certain loss function [24]. This problem was quickly solved [25], but Paul Davenport provided the first 
really useful solution [12]. He rewrote Wahba’s original loss function in quaternion form as 


• L = a t |b. - Ar t f =\- trace(4//) = \- q T Kq, 


(13a) 


where r ( is a vector to an observed object coordinatized in a reference frame, b ( . is the corresponding vector 
coordinatized in the spacecraft body frame , and <z. is the non-negative weight assigned to this observation; 
and where 

(13b) 

It follows that the optimal quaternion is the eigenvector of K with the maximum eigenvalue: 


\ = X; a t , B = 2/ a i b i rf , and K =■ 


B+ B T - / 3x3 trace B X, a t b f x r. 
(X,ab, xr.) 7 " trace .B | 


^opt ^'max^opt ‘ 


(14) 


This is Davenport’s q method, which is greatly appealing to physicists who love to work with eigenvectors 
and eigenvalues. Although very robust, Davenport’s q method is also very slow, and Shuster’s QUEST has 
been much more widely applied. 

I was aware of the Singular Value Decomposition (SVD) algorithm, but it was not until Malcolm Shuster 
applied it to a sensor calibration problem that it occurred to me to use this technique to solve Wahba’s 
Problem [26]. The SVD of the matrix B in Eq. (13b) is given as 

B = USV T , (15a) 

where S is diagonal and U and V are orthogonal. It was easy to show that the optimal attitude matrix is 

A opt = U diag([l 1 det U det V]) V T . (15b) 
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The SVD algorithm is as robust as Davenport’s q method, but even slower. It has provided useful analytical 
insights into the solution of Wahba’s problem, however. 

One drawback of all these solutions to Wahba’s Problem is that they estimate the attitude only, unlike more 
general estimators that can simultaneously solve for other parameters such as gyro drifts. I attempted to 
remedy this by considering the body and reference frame components of the observations to be functions of 
some parameters, and using a perturbation expansion in these parameters. The application of this theory did 
not live up to my hopes, but the papers contain some interesting algebraic identities that Ed Seidewitz and I 
discovered along the way [27, 28]. It was a great loss for the aerospace community when Ed left Goddard 
to pursue his first love of software engineering. 

Although the parameter estimation papers were not successful, they led to the discovery of a fast optimal 
attitude matrix (FOAM) solution to Wahba’s problem [29]: 

A opt = rV + HI*) 5 + K^ bT - BB T B ] , (16a) 

where • denotes the Frobenius norm, adj(-) denotes the classical adjoint, 

II IIF 

and >^-detS. (16b) 

As 1990 approached, Malcom Shuster proposed that we write a survey paper entitled “25 Years of Wahba’s 
Problem.” We could never find the time for this, but in 1995 I gave a lunchtime talk called, unsurprisingly, 
“30 Years of Wahba’s Problem.” This talk came to the attention of Daniele Mortari, who had invented 
several new solutions to the problem. He proposed that we write a survey presenting the different methods 
and comparing them for speed and accuracy. In preparing the paper, I developed first-order perturbation 
forms of Daniele’s ESOQ and ESOQ2 algorithms. We found that ESOQ and ESOQ2 are the fastest of 
currently known methods, although not significantly faster than the better-known QUEST algorithm [30]. 

Orbit Analysis 

About the time I went from NRL to Goddard, the radical idea was proposed that every engineer should 
have a computer on his or her own desk. Among the applications proposed to sell this idea to management 
was a simple mission analysis program to plot ground tracks of near-Earth satellites. This used a simple 
analytic orbit propagator incorporating the secular effects of the Earth’s oblateness. Comparison with a 
more precise propagator showed an in-track error increasing linearly with time. Remembering something 
Bob Dasenbrock had told me about semianalytic propagators, I gave him a call. He said that the errors 
arose from the difference between the mean and osculating semimajor axis, leading to an incorrect value 
for the mean motion. I added a linear mean-to-osculating transformation and its exact nonlinear inverse to 
the propagator, giving a zeroth-order extended-time algorithm. This fixed the problem, and the propagator 
was incorporated by Jim Jeletic into the Mission Planning Graphics Tool (MPGT) [31], and was also used 
to compute orbits for a wall display at the Smithsonian Air and Space Museum. 

An opportunity arose in 1989 to transfer to Henry Hoffman’s Guidance and Control Branch. Remembering 
my work on SMM as a high point of my career, I seized this opportunity to work with Hoffman, Donohue, 
and Flatley, and to get closer to the actual design of attitude control systems. 

Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX) 

The first spacecraft I worked on I after my transfer was SAMPEX, the first of the Small Explorer (SMEX) 
series of spacecraft. Its attitude control system had been largely defined by Tom Flatley before I entered the 
picture [32], It had a single reaction wheel providing an angular momentum bias along the spacecraft’s 
Sun-pointing y axis. Reaction wheel torques controlled the pointing of the z axis science instruments. 
Magnetic torquers were used to damp nutation and keep the y axis pointed at the Sun. The attitude was 
sensed by a three-axis magnetometer and a two-axis digital Sun sensor with a ± 64° field of view centered 
on the y axis and with 0.5° resolution. The TRIAD algorithm, or vector method [12], was used to find the 
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attitude matrix A. This required the spacecraft to carry onboard spacecraft and Sun ephemerides and a 
magnetic field model. The spacecraft ephemeris was a simple integration of F = ma , with a geopotential 
incorporating J 2 , J 3 , J 4 , S 2 2 , and C 22 corrections and a drag correction using a Jacchia- Roberts density 
model, reinitialized periodically with uplinked Cartesian orbit vectors. My innovation was adding S 22 , C 22 
and the Jacchia-Roberts model. The propagator exhibited in-track errors that were negligible for SAMPEX, 
but became problematic for more precise applications. Their cause and cure were not discovered until much 
later by Mark Beckman [33]. The cure is related to the solution of the errors in the MPGT propagator 
pointed out to me earlier by Bob Dasenbrock; the orbit must be put on the correct energy shell by 
initializing the propagation at a point where the approximate onboard gravitational potential is numerically 
close to the more exact potential used to compute the initial Cartesian orbit elements. 

The spacecraft angular velocity was found by 
differentiating the attitude matrix 

[a )X] = -AA T , (17) 

where [• x] denotes the cross-product matrix. Then 
the total system angular momentum was given by 

H = ^ + tf w heeli> ( 18 ) 

The magnetic control maintained the angular 
momentum with a magnitude of H 0 pointed 
simultaneously along the spacecraft y axis and at the 
Sun. An angular momentum error was computed as 

“ AH = (H - # 0 j) + (H - H 0 s) , (19) 

where s is the Sun unit vector in the body frame. This 
was driven toward zero by a magnetic dipole 

m = ^mag( AHxB )> (2°) 

where B is the magnetic field in the body frame. 



Figure 3. SAMPEX 


This simple control law worked well in simulations until the finite resolution of the digital Sun sensor was 
modeled. When this resolution was included, the TRIAD attitude solutions had discrete jumps, and the 
extremely noisy angular rates computed by Eq. (17) destabilized the attitude control. My solution was to 
add a simple constant-gain Kalman filter [34] 


H = 0 -^predicted + ^derived - 


(21a) 


where K is the filter gain, H deriveci is given by Eq. (18) and 

H predicted« = A(t)A r (t- At)H(t- Af) + (m x B)Af . (21b) 

SAMPEX was launched on July 3, 1992 with a lifetime requirement of two years and a goal of three years. 
It is still producing valuable scientific data, and its z axis pointing law has been reprogrammed several 
times to meet changing observational goals [35]. 

Geostationary Operational Environmental Satellite (GOES) 

In 1990 NASA and NOAA, the National Oceanic and Atmospheric Administration, began a study of an 
improved version of the Geostationary Operational Environmental satellite (GOES), which was itself an 
improvement on the SMS satellite that I had studied for Mel Velez many years earlier. Frank Bauer headed 
up the ACS portion of this GOES-Next study. My contribution was to develop an attitude determination 
concept to provide the required arcsecond pointing accuracy. The GOES-I/M spacecraft had pushed the 
accuracy attainable with Earth horizon sensors about as far as it seemed possible to go. The obvious (to me) 
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solution was a stellar-inertial system similar to that developed by Jim Murrell for MMS. The major effort 
was to find orientations for the star trackers that would protect them from Sun impingement and to use both 
Farrenkopfs analytic equations and numerical simulations to show that the required pointing accuracy 
could be attained. The idea of a Landsal ACS in geostationary orbit seemed obvious to me, but one 
reviewer for the Journal of Guidance, Control, and Dynamics , where the concept was published, found it 
very innovative [36]. It formed the basis for the winning proposal for the GOES-N/P spacecraft. 

Hubble Space Telescope (HST) 

HST was released by STS-31 in April 1990, less than a year after I joined Henry Hoffman’s branch. Before 
my transfer, I had worked with Paul Davenport on some of the optical calibration algorithms being 
developed for HST. This was the only time I worked with Paul, whose reluctance to publish prevented him 
from achieving the wider recognition he deserved. 

HST exhibited two problems almost immediately after launch. The most famous is the aberration of 
primary mirror, which was not a control problem. The second problem was the appearance of attitude 
disturbances that were most severe at entry and exit from the Earth’ s shadow, but that persisted throughout 
the daylight portion of the orbit. Henry immediately attributed these disturbances to the bistem booms used 
in the solar arrays, which had a tendency to twist as shown in Figure 5. Henry, Tom Flatley, and Jim 
Donohue immediately set to work to develop a disturbance-attenuating control algorithm for HST. 
Engineers at Marshall Space Flight Center and at Lockheed-Martin, the HST prime contractor, also 
attacked the problem. This work led to the Solar Array Gain Augmentation (SAGA) algorithm that HST 
employed until the solar arrays were replaced with rigid arrays in 2002. 1 was largely an admiring spectator 
of this development, but I did play a role in the development of a second- generation SAGA-II algorithm. 


Figure 4. GOES-N 


Figure 5. Hubble Space Telescope 
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The next problem to surface on HST was the failure of two of its six precise air-bearing gyros. Since HST 
needs three gyros for control, failure of two more gyros would have actuated an undesirable firmware 
controller using limited-life ball-bearing gyros, so the HST Project decided in 1991 to develop a software 
safe pointing mode requiring fewer than three gyros. As Goddard’s lead engineer for this effort, I worked 
with a contractor team headed by John Nelson of Lockheed-Martin. Their preliminary study had shown that 
a zero-gyro safe pointing mode could be implemented, and I quickly agreed that this was preferable to a 
two-gyro or one-gyro mode. We developed the Zero-Gyro Sunpoint (ZGSP) mode from initial concept to 
flight readiness in four months [37]. This mode points a preferred axis toward the Sun during the daylight 
part of the orbit and uses a momentum bias in the reaction wheels along this axis to hold attitude during 
eclipse. This is similar to the science mode of SAMPEX, but ZGSP does not rely on an onboard ephemeris 
or magnetic field model, so TRIAD is out of the question. Instead, the ZGSP mode controls rotation around 
the Sun line by locking it to the Earth’s magnetic field. This has the unfortunate property that for some 
orientations of the orbit plane, magnetic equator, and Sun line, the apparent rotation of the inertial magnetic 
field as HST orbits can cause the body angular momentum to cancel a significant fraction of the bias 
momentum. The ZGSP mode has exhibited some large attitude excursions arising from this effect during 

eclipse, but it maintained HST in a power and thermally safe state for 38 consecutive days in November _ 

and December of 1999 after the failure of four gyros and until the replacement of all six gyros. 

It warmed my physics- trained heart that Henry Hoffman always emphasized the importance of the 
conservation of angular momentum. He was quick to discredit any analysis that reached, erroneous 
conclusions by ignoring this principle, something that I once experienced personally to my great 
embarrassment. I was delighted to find an application of angular momentum conservation to onboard 
detection of a spacecraft failure. The HST computer calculates the total system angular momentum as the 
vector sum of reaction wheel momentum sensed by the wheel tachometers, and spacecraft body momentum 
sensed by the gyros. The high torque capability of the reaction wheels can cause both of these components 
to change rapidly, but the total system momentum only changes slowly as a result of environmental 
disturbance torques. Propulsion torques would violate this assumption, but HST has ho propulsion system. 
Thus a rapid change in system momentum would indicate a failure of either a reaction wheel tachometer or 
of a gyro. A reaction wheel tachometer failure can be identified on HST by an independent test, a 
comparison of the reaction wheel momentum change to applied torque, so the system momentum test was 
implemented to identify gyro failures [38]. It has not identified any failures to date, however. 

Research in Attitude Estimation 

In his last year at Goddard, Gene Lefferts had a National Research Council Resident Research Associate 
(NRC RRA), S. Vathsal, whom I inherited after I replaced Gene. Vathsal developed a second-order Kalman 
filter based on the methods presented in the Lefferts, Markley, and Shuster paper [39]. I have had the 
privilege and the pleasure of working with several other NRC RRAs, including Itzhack Bar-Itzhack from 
1987 to 1989, Nadav Berman from 1992 to 1994, and Yaakov Oshman from 1996 to 1998. Although I was 
nominally their advisor, I learned more from them than they from me, and our collaborations produced 
some interesting research that I certainly could not have accomplished without them [40-42], 

Working with these NRC RRAs, most especially with Yaakov Oshman, increased my understanding of 
spacecraft attitude estimation — deepening my knowledge of estimation theory, opening my eyes to a 
variety of approaches, and challenging my assumptions and approximations. When I finally thought I 
understood the issues well enough, I wrote an update the Lefferts, Markley, and Shuster paper [43]. 
Extensive discussions with Russell Carpenter at Goddard sharpened my arguments and significantly 
improved the clarity of the final paper. Russell also guided me into the mysteries of Daum’s approach to 
estimation [44], which contributed- greatly to my most recent work on attitude filtering [45]. 

I have saved my most rewarding NRC RRA experience until last John Junkins’ student Joe Mook had 
presented some of their work on Minimum Model Error Estimation at a Goddard Flight Mechanics/ 
Estimation Theory Symposium in 1988 [46], All of John Junkins’ many graduate students that I ran into 
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did interesting work, so I began to regard this as a theorem. I saw Joe periodically at conferences, and we 
became friends. His reward to me was to send his best graduate student, John Crassidis, to be my NRC 
RRA for 1994-96. 1 quickly discovered how creative and prolific John is; during his stay at Goddard he 
presented more than a dozen conference papers and submitted a half-dozen for publication, with me as a 
co-author. An example of my contribution was to sketch a method for approximate GPS attitude estimation 
that John worked out in a week while I was away, along with an alternative method, a complete covariance 
analysis, and'a comparison with the covariance of the optimal estimate [47]. Our collaboration did not end 
with the termination of John’s tenure at Goddard, but continues to this day, much to my delight 

I also returned to Farrenkopf’s analysis of the steady-state single-axis attitude estimation errors using gyros 
and angle measurements. Spacecraft generally use rate-integrating gyros, which have random error on the 
angle output with standard deviation & e that he did not include [12]. I discovered that including this term 

led to a quartic equation that could be solved exactly like Farrenkopf’s. Using the same notation as in Eqs. 
(10)— (12) and defining 

S e= a e/°n’ S u 3 ( At f P ° u! ° n > 311(1 ^ = ( At) 1 /2 £T y /<T„ , (22) 

the steady-state errors are 

W = -C^(T„(A tf 12 , p ge (~) = (£ 2 - 1 )a 2 n , P eg (+) = (1 - r V„ 2 , 

and P hh (+) = cr„K 2 +2y a u a n (At) 1 /2 + 1 - a 2 u (Atf ] l/2 ± \<J 2 U At , (23) 

where C = 7 + K + ±(2 r S a + S 2 + I,S 2 ) l/2 and y = (l + S 2 +±S 2 + ^S„ 2 )- I £ ; (24) 

These results are identical to Farrenkopf’s if (7, = 0 . 1 showed my original messy derivation of these 

equations to Reid Reynolds, and together we found the much more straightforward path to the results that . 
appears in the published version [48], 

Tropical Rainfall Measuring Mission (TRMM) 

In 1991, GSFC embarked on the simultaneous development of two spacecraft in-house: the Tropical 
Rainfall Measuring Mission (TRMM) and the X-Ray Timing Explorer (XTE, later renamed RXTE in honor 
of Bruno Rossi). The two spacecraft used a common procurement for reaction wheels and gyros, but their 
attitude control systems were quite different. XTE chose a stellar-inertial reference designed by Mike 
Femiano based on Murrell’s MACS design to satisfy its arcsecond pointing requirements [49], while 
TRMM chose an Earth-referenced system to satisfy its less stringent 0.3° (3cr) attitude knowledge 
requirement. TRMM employed a static Earth horizon sensor similar to that used on the NOAA low-Earth- 
orbiting spacecraft and on the Defense Meteorological Satellite Program (DMSP) spacecraft. In early 1994, 
after development of the TRMM ACS was well underway, we became aware of a progressive loss of 
optical transmission through of some of the windows of the horizon sensors on the DMSP spacecraft. We 
worried that TRMM, with its lower 350 km altitude orbit, might experience a more severe problem, and 
possibly even complete loss of the attitude reference. To protect against this possibility, we considered 
adding a star tracker to TRMM, which would have incurred severe monetary and schedule penalties. 

Recalling that Joe Hashmall at CSC had found some promising results on batch attitude determination 
using magnetometers and gyros, I sought more information from Eleanor Ketchum of the Flight Dynamics 
Analysis Branch, who was sponsoring his work. They and Joe Sedlak responded enthusiastically, and 
within a month had dev eloped.a Kalman filter that gave 0.3° attitude accuracy using real data from the. 
Upper Atmosphere Research Satellite (UARS) [50]. The common procurement of XTE and TRMM gyros 
permitted this filter to be used on TRMM, since the gyros were much more accurate than required by the 
Earth-sensor-based TRMM ACS but provided the dynamics memory of approximately one orbit period 
required by the filter to reduce the effects of magnetic field modeling errors. 
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John Crassidis, Kong Ha, and I studied alternative 
Kalman filter implementations using magnetometer, 
Sun sensor, and gyro data [51]. Although various 
simpler algorithms were promising, the TRMM 
ACS group decided to adapt the XTE Kalman filter 
for this purpose; and this was accomplished under 
the direction of Steve Andrews, In 2001, it was 
realized that TRMM had insufficient propellant to 
continue compensating for atmospheric drag in its 
350 km orbit, so the orbit was raised to 402 km. The 
static Earth sensor that had performed according to 
specification at the 350 km altitude ceased providing 
valid data above 380 km, so the Kalman filter was 
enabled. After some retuning, the filter is providing 
attitude accuracies of approximately 0.2°, about the 
same as the horizon sensor had provided [52], 

Wilkinson Microwave Anisotropy Probe (WMAP) 


Figure 6. TRMM 


The MAP proposal effort began in late 1994, anticipating a Mid-sized Explorer (MID EX) Announcement 
of Opportunity in early 1995. The Principal Investigator (PI) was Chuck Bennett of Goddard, who had been 
deputy PI on the Differential Microwave Radiometer (DMR) instrument on COBE. The MAP proposal was 
a Goddard-Princeton collaboration to repeat the DMR measurement of the cosmic microwave background 
(CMB) with increased resolution and precision. To that end, MAP was to be launched to the Earth-Sun L2 
libration point to avoid the Earth’s atmosphere, radiation, and magnetic field. The Princeton team included 
David Wilkinson, the dean of CMB studies, in whose honor MAP was rechristened WMAP in 2002. It was 
a great privilege to participate in planning meetings with him. Cliff Jackson, who headed the Goddard 
engineering team and later became the MAP Instrument Systems Engineer, deserves a lot of credit for the 
success of the proposal. 

The scientists wanted to obtain a highly interconnected set of measurements over an annulus between 45° 
and 90° from the anti-Sun line, requiring MAP to execute a fast spin at about 0.5 rpm and a slow precession 
of its spin axis at one revolution per hour at a constant angle of 22.5° from the Sun line. I realized that this 
scan pattern could be easily represented by 3-1-3 Euler angles and could be accomplished by a zero- 
momentum ACS controlled by reaction wheels, just as for COBE. John Crassidis took the Euler angle 
description and used MATLAB to produce the “spirograph” shown in Figure 7. The science team were 
delighted that we could produce this pattern; and a modified version of Figure 7 was included in the 
proposal. Tom Flatley tried without success to find a dual-spin dynamical configuration that could provide 
this scan pattern; and I believe that if he couldn’t, nobody could. MIDEXs were to be single-string 
spacecraft, so the MAP ACS concept had only three reaction wheels to serve the dual function of 
counterbalancing the body’s spin angular momentum to maintain zero system momentum and applying 
control torques to provide the desired attitude. The wheel axis orientations were chosen to bias all wheel 
speeds away from zero to avoid undesirable zero-speed crossings. 

Attitude sensing was by a star tracker and a digital Sun sensor, despite the preference of another Goddard 
engineer (not in Hoffman’s branch) to use an Earth sensor in place of a star tracker. My GOES and TRMM 
experience had left me with a prejudice against Earth sensors, and the Earth/Sun separation of less than 10° 
as viewed from the MAP orbit at L2 would have led to marginal attitude determination accuracy. I was 
relieved to find two star tracker manufacturers who could provide a tracker to acquire and track stars at 3 
deg/sec, as required by MAP. Ironically, the star tracker contract was eventually awarded to a third vendor. 
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Figure 7. MAP Scan Pattern Figure 8 . MAP 


I had laid out this basic ACS concept before the end of 1994, and John Crassidis developed a simple 
MATLAB simulation to confirm its viability [53]. MAP was selected in 1996 to be one of the first two 
MIDEX missions. The ACS was developed by a MIDEX ACS group already in place at Goddard, headed 
by David Ward. A two-wheel backup mode was later developed to provide reduced performance in the 
unlikely event of a reaction wheel failure, which has not occurred [54]. 

MAP was launched in 2001, and its data have enabled scientists to determine the values of key 
cosmological parameters and to answer questions about the origin of structure in the early universe and the 
fate of the universe. 

SUMMARY 

My career did not end with MAP’s launch in 2001, although I still see MAP as its high point. I could 
mention Goddard’s Swift spacecraft, a MIDEX mission launched in 2004, that uses six reaction wheels to 
slew 60° in 60 sec to catch elusive short-lived gamma-ray bursts on the fly. I continue to follow HST, 
which is currently observing in two-gyro fine-pointing mode, saving the remaining two operable gyros until 
a new servicing mission can replace all six gyros, and with a one-gyro science mode under development. I 
am anxiously awaiting the launch of GOES-N, held up since last winter by a variety of difficulties, most 
recently a labor dispute. I have largely been an interested observer and friendly critic of these 
developments, however. The torch has been passed to a new generation. 

My aerospace career has been extremely satisfying, but I will especially cherish the Dirk Brouwer Award 
because of the high regard in which I hold you, my colleagues in the AAS. Your opinion means the world 
to me. I look back on my career as having been largely fortuitous; a matter of finding myself in the right 
place at the right time. It has of course also required me to make the best use of my opportunities when I 
found them. 

My greatest reward has been the opportunity, it has given me to work with incredibly talented colleagues. .... 
Many of the people named in this paper are well known to you, including three previous recipients of the 
Dirk Brouwer Award, while others have achieved only more limited reputations. I hope that my talk will 
serve to rescue some of these valued colleagues from undeserved obscurity. 
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